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Abstract 

We investigate the quantum fluctuations of a single atom in a weakly driven cavity, where the 
center of mass motion of the atom is quantized in one dimension. We present analytic results for 
the second order intensity correlation function g^'^^r) and the intensity-field correlation function 
hg^r), for both transmitted and fluorescent light for weak driving fields. We find that the coupling 
of the center of mass motion to the intracavity field mode can be deleterious to nonclassical effects 
in photon statistics; less so for the intensity- field correlations. 
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INTRODUCTION 



Since the mid 1970's, quantum opticians have been investigating exphcitly nonclassical 
states of the electromagnetic field, and ways to determine if a field state is nonclassical. 
These types of states are ones for which there is no underlying non-singular probability 
distribution of amplitude and phase, or more technically, they exhibit a positive definite 
Glauber-Sudarshan P distribution. Much work has focused on photon antibunching, sub- 
Poissonian photon statistics, quadrature squeezing, and entangled atom-field states The 
generation of such light fields may have applications in quantum information processing, 
atomic clocks, and fundamental tests of quantum mechanics, for example. One system that 
has long been a paradigm of the quantum optics community is a single-atom coupled to a 
single mode of the electromagnetic field, the Jaynes-Cummings model|2]. In practice the 
creation of a preferred field mode is accomplished by the use of an optical resonator. This 
resonator generally has losses associated with it, and the atom is coupled to vacuum modes 
out the side of the cavity leading to spontaneous emission. Energy is put into the system 
by a driving field incident on one of the end mirrors. The investigation of such a system 
defines the subfield of cavity quantum electrodynamics [3]. The presence of the cavity can 
also be used to enhance or reduce the atomic spontaneous emission rate [3]. This system has 
also been studied extensively in the laboratory, but several practical problems arise. [H O |6] 
There are typically many atoms in the cavity at any instant in time, but methods have 
been developed to load a cavity with a single atom. A major problem in experimental 
cavity QED stems from the fact that the atom(s)are not stationary as is often assumed by 
theorists. The atoms have typically been in an atomic beam originating from an oven, or 
perhaps released from a magneto-optical trap. This results in inhomogeneous broadening 
of the atomic resonance from Doppler and/or transit-time broadening. Using slow atoms 
can reduce these effects, but the coupling of the atom to the light field in the cavity is 
spatially dependent, and as the atoms are in motion, the coupling is then time dependent; 
also different atoms see different coupling strengths. 

With greater control in recent years of the center of mass motion of atoms, developed by 
the cooling and trapping community, preliminary attempts have been made to investigate 
atoms trapped inside the optical cavity [Tj. The recent demonstration of a single atom laser 
is indicative of the state of the art [H]. In this paper we consider a single atom cavity QED 
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FIG. 1: Single atom in a weakly driven optical cavity with an external potential 



system with the addition of an external potential, provided perhaps by an optical lattice, 
and study the photon statistics and conditioned field measurements of both the transmitted 
and fluorescent fields. We seek to understand (with a simple model at first) how the coupling 
of the atom's center of mass motion to the light field affects the nonclassical effects predicted 
and observed for a stationary atom. 

The system we consider is shown schematically in Fig. 1. We utilize the quantum tra- 
jectory method in which the system is characterized by a wave function and non-Hermitian 
Hamiltonian 
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where we also have collapse operators 
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associated with photons exiting the output mirror and spontaneous emission out the side 
of the cavity. The indices e{g) indicate the atom in the excited (ground) state, n is the 
photon number, and / is a quantum number associated with the presence of bound states 
of the external potential. We have the usual creation (a) and annihilation (a"!") operators 
for the field, and Pauli raising and lowering operators a± for the atom. The bare energies 
are -En,z,g = fi{nuj + fli) and -En,z,e = fi{{n + l)uj + fli), where the fli are the discrete, 



bound, states of the external potential. The classical driving field (in units of photon fiux) 
is given by Y. We take the external potential in which the atom is trapped to be harmonic 
along the cavity axis, Vext = cn{z — zo)'^/2, which could be appropriate for a 1-D optical 
lattice inside the cavity. We ignore the generally weak transverse dependence of the atom- 
field coupling, g{f) g{z) = gmf {k{z — zq)), with the maximum coupling given by g^. = 



IJieg\Juj /2htQVaxid f {k{z — zo))is the cavity field mode function. Here /leg is the dipole 
transition matrix element, and V is the volume of the cavity mode. We assume for simplicity 
that the bottom of one of the lattice wells coincides with an antinode of the cavity field. 
Following the treatment of Kimble and Vernooy. PH] and keeping only terms to {z — Zq)"^, we 
find the non-dissipative parts of the Hamiltonian to be 
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H = ^ + -{z- zo^ + inE{a^ - a) 
2m 2 



+hgm{l + ^^^^] (aV_+aa+) (5) 



where the characteristic distance rj is defined by 



nd?giz) _ d?fiz) 

For a standing wave mode, f{k{z — zq)) = cos{k{z — zq)) we have 1] = = X/27T, where A is 
the wavelength of the cavity mode. Consider the action of this Hamiltonian on the dressed 
states I n, ±) = (l/y^) (| n,g)± \ n — l,e)) with n the number of intracavity photons, and 
e{g) denotes the excited (ground) state of the atom. As (aV__ + ao"+) | ra, ±, /) = ±-\/n | 
n, ±, I) we have an effective potential 



V{z) = \L± j {z - z,r = Imnl^ (7) 



with an effective harmonic frequency Qn± defined above. In the dressed state basis, the 
selection rule for dipole transitions isAZ = 0. It is worth noting that in a basis defined by 
an outer product of the atom-field dressed states and the bare vibronic levels of the external 
potential enumerated by the quantum number L, which we call the casually dressed states, 
the selection rule is AL = 0, ±2. These arise from absorption of a photon traveling to the 
right (left) in the cavity, with reemission into the same direction (AL = 0), while absorption 
of a photon traveling in one direction and emission into the opposite direction leads to a 
momentum kick for the atom or 2hL, leading to AL = ±2. 
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FIG. 2: Energy level diagram 

We then use the dressed states | n, I, ±) where the index I denotes the degree of excitation 
in the vibronic states corresponding to combined lattice/field couphng potential. Please note 
that these are not the vibronic states of the optical lattice alone. An energy level diagram 
is shown in Fig. 2. We notice that the level spacings of the three sets of dressed vibronic 
states are not equal, due to the ±gy/n term in the vibronic frequency. 

INTENSITY-INTENSITY CORRELATIONS 

We next consider the second-order intensity correlation function (7*-^-'(t) = 
(at(0)at(r)a(r)a(0))/(ata)|^. 

In the weak field limit, only states with 2 or fewer quanta of energy are left within the basis 
(we must keep states with at least two photons, as we wish examine photon coincidences). 
The limit we are considering is one in which Y ^ 0; and we truncate our equations of 
motion to lowest order in Y. If no driving field were applied, the atom would certainly be 
in the ground state, so we make the approximation that for weak fields Co,i.g ~ 1- With 
no trapping potential, one would have Cg^ = 1.0; here we must specify the set of initial 
populations which correspond to the center of mass motion of the atom, subject to the 
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normalization condition J2i |C'o,«,gP = 1-0. The potential is taken to be of the same sign 
for plus and minus dressed states, which is possible by placing the lattice field at a "magic 
"frequency [HI [T21 [13] -The driving field is responsible for populating the atom's excited 



states, and thus Ci,/_± ~ YCo^i^g ~ Y. 
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reasoning 


; can be continued and we determine 


that our scaling should be 
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In the weak field limit, the one excitation amplitudes satisfy 

'7 K \ Y 
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Ci,i,- (9) 

Ci,i,- = - ( 4 + 2 + '^'-i-) ^i'''- " 

^ Ci,i,+ (10) 
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with Ai^i± = (f2i,± — i^o) ^ =t (7, recall the effective harmonic frequency 

and again, we keep only lowest order terms in the driving field Y . These are the frequencies 
that correspond to the energy levels of the system, the a term arises from the external 
potential, the -^g^Jn terms from the spatial structure of the cavity mode function and the 
coupling of motion in the mode to the interaction with the driving field. 

As a first step we note that we can solve the equations of motion for the slowly varying 
population amplitudes Dn,k,±,, defined as 
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We find that our Cn,i,± equations become 

Do,i,g = Co,i,g - 



(13) 
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In the weak field limit, these equations have a steady state solution 

= E fe+l^' +> + Dn.i,-\n, I, -)) 



(15) 



(16) 
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The system reaches this steady state, which has a very small average photon number, 
(a^a) ~ <^ 1. In any given time step At the probability of a collapse is given by 



Pcav = 2n{o) a) ^ 1. Similarly the probability of a spontaneous emission event in a time 
step At, Patom = l{<^+o'-)^t is small. Eventually there is a cavity emission, or a spontaneous 
emission by the atom, leaving the system in the states 




I^ct(O)) 



Transmission 



" ' = ||l^^jff|2 : Flourescence 

In the steady state, all population amplitudes Dn,k,± are constant, and we may set all 
Dn^k,± = 0. Equations 14 and 15 then become 
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Solving for D^^^ ^ and D J'^^ _ , we find 



where 
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FiHi-G^ ' 



(21) 



^ + ^ + - Qoi - g) 
^ + + i{Qi^+l ~ QqI + g) 



Using the same procedure, we may use our results for Df^i _,_ and Df'^i _ and solve 

—G(32i — ZiPii 



(22) 

equations 



(16) and (17) for Dfi^ and -D|*;_, finding 
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where 



Yi = ^ + y + z(n2,+/-^^o^ + v^^7), 
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P2,l = +Y 
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Now that the steady state values for the population amplitudes have been calculated, our 

task is to solve for the time evolution of ; and Di^i _{t). The probability of a cavity 
emission at time r given that one occurred at r = 0.0 is 2K{'^cTWo'\^CT)^'t-i hence we have 
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(25) 



where we have truncated the results to lowest order in the weak field limit. Similarly we 
have for the second order intensity correlation function for the fluorescent fleld is given by 



{^ss\<y+<y-\^ss) 

(26) 

To facilitate solving the time evolution of the one-excitation amplitudes we write them in 
matrix form as 

A{t) = MA{t) + A, (27) 
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where 



A{t) = 
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The form of the time evolution of and 

A{t) = (^e^*^-^) A(0) + (Se^'S-^ - l) M'^A. 
Without showing the details of such calculations, we arrive at 



(29) 
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Xi = 



X2 



Fi + Hi 



^{Fi-HiY + iG^ 



Ai = X1 + X2, 

A2 = Xi - X2, 

^2^ Fi- A2, 

= FiHi-G'. 



(32) 



We are now equipped with all the necessary information to solve for the dynamics of our 
system in the weak field limit. 

For the rest of the paper, we restrict ourselves to the deep trapping limit, where a > 
gmy/n/X^. We may then use the binomial approximation, and define 



fal hmgmVn 
V m Zrj'^a 



(33) 
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therefore we find A„^± = Qn, ± — for n = 1, 2 to be 

Ai,_ = -Ai,+, 
A2,+ = V2Ai,+, 

A2,_ = -V2Ai,+. (34) 

and we can characterize everything by the one detuning Ai__|_. This is analogous to the 
Lamb-Dicke regime in an ion trap. 

By using the well dressed states, we have a set of equations that will have a steady-state. 
Recall that the quantum number I is associated with a well-dressed state, and not simply the 
vibrational quantum number of the lattice potential. To solve these equations it is necessary 
to specify the amplitudes, Cc,,o.i(0) that are each of order unity. They can be related to the 
initial center of mass state of the atom via C(,,o,z = (5',0,/|'^) = / (j)*gQi^{x)dx , or simply 
specified. 

For weak driving fields, the probability of a cavity emission in a time Ai is given 
by Pcav — 2K{a^a)/S.t is quite small, as is the probability of a spontaneous emission, 

Pspon.em. 

= 7(o-+o-_)At. In this case the wave function attains a steady state 1^)55 = 

00 r -1 

E Cff+e-*^i-'>+*|l,/,+) + Cff_e-*^i.'.-*|l,/,-) . 

n,l=0 L ' ' ' ' J 

After a photon is detected in transmission, at t = the wave function collapses to 

mO))coii = a\^)ss/\a\^)ss\ (35) 
00 

n,l=0 

+C3(t)e-^--'*|e,n,/)) (36) 

. The initial value of the one-photon amplitudes of the collapsed state are related to the 

steady state two-photon amplitudes 

^g,iA^) - ^ [91^55 \^,\rss PI ^^'^ 

l^n,l [^1 <^g,2,i I + I <-^e,l,/ I 

Cfol^(O) = ^ (38) 

l^n,l [^\ ^g,2,l I + I ^e,l,l 

and these are found above. 
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The relation between the well dressed probability amplitudes and the bare amplitudes is 

Dn,l,± = -^{Ce,n-l,l ± Cg,n,l) (39) 

Before turning to our results, let us recall the relations that g^'^^r) must satisfy if the 
field can be described by a classical stochastic process; if it has a positive definite Glauber- 
Sudarshan P distribution, 



9^'\0) > 1 (40) 
^(2)(o+) > ^{2)(o) (41) 

|^7^'^(r)-l| < \9^'\r)-l\ (42) 
Violations of all three of these inequalities has been observed in CQED systems [5l [6] 

RESULTS FOR INTENSITY CORRELATIONS 



As the system has a steady-state wave function in the steady state, we may write 

r)\a^a 

where we define 



^(2)(^) ^ {^cir)\a^a\^ciT)) _ ^^^^ 



1^^) = "'^^^^ = ^ (44) 

l«l^ss)| ^{\'^ss\a^a\^ss)\ 

In Fig. 3 we plot 9^'^\0) — 1 for an initial state | where there is equal population in the 
I 0, /, g) states for I = 0, Imax- As more states are involved, we find that the antibunching goes 
away. This is due to the fact that the two single-photon vibronic ladders have a different 
frequency spacing than the ground state vibronic levels, and is consistent with the effect of 
detunings on the photon statistics Involving more / states makes the width of the state 
larger, increasing Ax for the center of mass wave function of the atom. The antibunching 
also goes away if we just prepare the system in a particular higher / state | 0,lo,9)- The 
optimum state would seem to be the ground state of the bare vibronic potential. 

Instead of just assigning values to the probability amplitudes (subject to normalization) 
we can specify the center of mass wave function and calculate the probability amplitudes 
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FIG. 3: g^'^\0) — 1 vs. Nmax, the highest occupied phonon number 

via Dg o i{0) = (\t'c'A/|0- If choose the center of mass wave function to be a Gaussian of 
width a, we can calculate these amplitudes easily using 



Dg,Q,l = {^Cm\1) 

= A 



oo 



(45) 



where y = x/a^ and ctq = y^/'^^o,i is the width of the ground state of the vibronic potential, 
and the normalization constant is A^^i = {m^lo^i/nhy^'^ /^/2"-n\. In Fig. 4 we show a plot of 
g^'^\0) as a function of cr/ao for parameters for which there is nearly perfect antibunching 
in the absence of an external potential. We see that there is a relatively wide region where 
the antibunching persists, but for cr/ao less than 0.2 or larger than 4, the antibunching 
vanishes completely. This can be understood by considering that a Gaussian wave function 
is superposition of various vibronic states, and that population of higher excited vibronic 
states is deleterious to antibnching. Only when a/cio ~ 1 do we have a center of mass wave 
function that has population predominantly in the ground state In Fig. 5, we show a plot 
of g^'^^r) for g/'j = 2, k/'j = 5, Ai +/7 = 0.1. Fig. 5a is for the atom initially in the ground 
state of the potential. We see that g^'^\0) is about 1.1. Classically, g^'^^r) could not then go 
below 0.9, but here it goes to zero. We refer to this as an undershoot. In Fig. 5b, we exhibit 
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FIG. 4: §^"^^0) — 1 vs. cr/co, the relative width of a Gaussian center of mass wave function 

g^'^^r) for an equal admixture of the ground state and fifth excited state. Here g^'^\0) is 4, 
and hence the fact that g^'^^r) is later zero is not nonclassical. The physical reason for this 
can be traced back to the fact that an atom in an excited state of the external potential 
is essentially detuned from resonance. With Ai_+/7 = 0.1, the detuning A^^^/j = 0.5. 
Previous work has shown that usually a detuning of half a linewidth is quite deleterious to 
nonclassical effects in g^'^^r). In Fig. 5c, we have results for what we refer to as a pseudo- 
Boltzmann. Here we populate 20 vibronic levels of the external potential at a "temperature" 
of 3mK. There is no decoherence associated with this distribution, i.e. all the off-diagonal 
matrix elements are not zero. This essentially results in a distribution over populations with 
small population in the first excited state, even less in the second, and so on. Here we see 
that with most of the population in the ground state, we essentially have the ground state 
result. In Fig. 5d, we show g^^^r) for an equal population in all 20 vibronic states. Here 
we see large photon bunching, and no nonclassical effects at all. This can be understood in 
terms of detunings of the various atomic states; this type of distribution over vibronic states 
would correspond to an atom more localized than the ground state of the external potential. 
Hence we see that localizing the atom too much results in a large spread in momentum 
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FIG. 5: Plots are g^^^ir) for 5/7 = 2, n/^f = 5, Ai,+/7 = 0.1, and (a) |0) only, (b) ^ [|0) + |5)], 
(c) Pseudo-Boltzmann, and (d) for 20 states with equal population. 



states that destroys the nonclassical effects. In Fig. 6, we see that the fluorescent intensity 
correlations are relatively insensitive to the choice of atomic center of mass wave function, in 
that ^''•^''(O) is due to the nature of single-atom fluorescence. In Fig. 7 we examine g^'^^r) 
for parameters where the transmitted intensity correlation function g^^\0) = 0.0. We see 
that for a superposition of ground and fifth excited states, we still have nonclassical effects, 
as g^'^\0) < 1. The initial slope of g^'^^r) is negative though, which is not nonclassical. 
For the pseudo-Boltzmann distribution, we see both types of nonclassical behaviors. In the 
case of equal population over 20 vibronic states, there is no nonclassical behavior at all. In 
Fig. 8, we look at a case where there is strong coupling, but no nonclassical behavior in 
the ground state case. We do have strong vacuum-Rabi oscillations. For an admixture of 
states, we see a beat frequency in the oscillations. The pseudo-Boltzmann case again is very 
similar to the ground state case. The oscillations are almost completely washed out when 
we have equal population in 20 vibronic states. In Fig. 9 we again look at a situation where 
the ground state case shows strong vacuum-Rabi oscillations as well as all three nonclassical 
behaviors; g^'^\0) < 1, the initial slope is positive, and there is an overshoot violation. The 



latter refers to g^^^'^r) violating the upper limit of the inequality in 42 In Fig. 10 we have a 
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FIG. 6: Plots are g^p^r) for 5/7 = 2, K/7 = 5, Ai,+/7 = 0.1, and (a) |0) only, (b) ^ [|0) + |5)], 
(c) Pseudo-Boltzmann, and (d) for 20 states with equal population. 
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FIG. 7: Plots are g^^ri'^) ^/t = 2-2, k/7 = 10, Ai,+/7 = 0.1, and (a) |0) only, (b) ^ [|0) + |5)], 
(c) Pseudo-Boltzmann, and (d) for 20 states with equal population. 
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) 

FIG. 9: Plots are gg^Cr) for 5/7 = 3, K/7 = .1, Ai,+/7 = 0.1, and (a) |0) only, (b) ^ [|0) + |5)], 
(c) Pseudo-Boltzmann, and (d) for 20 states with equal population. 
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FIG. 10: Plots are g!^liT) for 5/7 = 1, k/j = .77, Ai,+/7 = 0.1, and (a) |0) only, (b) ^ [|0) + |5)], 
(c) Pseudo-Boltzmann, and (d) for 20 states with equal population. 



situation where we only have an overshoot violation in the ground state case. This violation 
vanishes in the case of a superposition of ground and fifth excited state, as well as for an 
equal population of 20 vibronic states. In Figure 11, we examine the effects of increasing 
spacing between the vibronic levels. To this point we have dealt with detunings on the 
order of 0.1 linewidths. In Fig. 12 we can see that increasing the detunings allows us to 
see a larger effect due to the beat frequency. Changing the detuning to 0.3 and 0.5 of 7, 
we see that the initial slope is not nonclassical, but we still have ^'^^•'(O) < 1, and there is 
an undershoot violation. So the nature of the nonclassicality is not changed. At a detuning 
of 2.0, we still have an undershoot violation as well as evidence of oscillations at the beat 
frequency. In Figure 12, we examine the effects of increasing spacing between the vibronic 



levels. In this case we have antibunching, a violation of inequality in Eq.(40). Changing 



the Ai^+/7 to 0.3 and 0.5, we see that the initial slope is not nonclassical, but we still have 
(y((^)(0) < 1, and there is an undershoot violation[6j. So the nature of that nonclassicality is 
not changed. At a detuning of 2.0, (Fig. 5c) we still have an undershoot violation as well as 
evidence of oscillations at the beat frequency. 
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FIG. 11: Plots are grpJp{T) for g/^ = 1, k/^ = 1.6. All trials use [|0) + |5)] as the vibrational 
state distribution with (a) Ai,+/7 = 0.1, (b) Ai,+/7 = 0.3, (c) Ai,+/7 = 2.0, (d) Ai,+/7 = 0.5. 

WAVE-PARTICLE CORRELATIONS 



Recently, Carmichael his co-workers have introduced a new intensity-field correlation 
function /^^(t) that is of great interest [TH|15]. Because is an intensity-field correlation 
function, it takes the general form 

and for a quantized field, this becomes 

(at(0)a,(r)a(0)) 
^'^"^ = (ata)(ao) ' ^^^^ 

where we have, like for g^'^^r), exploited normal and time ordering, and we have used the 
quantum mechanical field quadrature operator : 



ae 2 



1 (^ae-'S ^ ^tg^^^) . (48) 



In Eq. (48), 9 is the phase of the local oscillator (LO) with respect to the average signal 
field. We see that with the a acting to the right, and the acting to the left at t = 0, a 
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FIG. 12: Plots are g^pj'ir) for g/j = 1, k/^ = 1.6. All trials use [|0) + |5)] as the vibrational 
state distribution with (a) Ai,+/7 = 0.1, (b) Ai,+/7 = 0.3, (c) Ai,+/7 = 2.0, (d) Ai,+/7 = 0.5. 



collapsed state is prepared, the collapse being a photon loss from the field, corresponding to 
a detection event. Then at t = r one measures (ag) conditioned on the previous detection. 
This differs from a direct measurement of (ag) with no conditioning. An ensemble average 
of the latter measurements (necessary to get a good signal to noise ratio) would yield zero 
due to phase fluctuations. The conditioned BHD measurement essentially looks at members 
of the ensemble with the same phase, a phase that is set by the photodetection. 

As with other correlation functions, like the second-order intensity correlation function 
g^'^^r), restrictions can be placed on hg^r) if there is an underlying positive definite prob- 
ability distribution function for amplitude and phase of the electric field, i.e. that the field 
is classical albeit stochastic. By ignoring third-order moments (a Gaussian fluctuation as- 
sumption that is valid for weak fields), one finds 

(: Aag{0)Aag{T) :) 



hg(T) = 1+2- 



(49) 



(AatAa) 

and we see that the intensity-field correlation function is connected to the spectrum of 
squeezing [H] 

Sg{uj) OC / dTCOs{uJT) [hg{T) — 1] . (50) 

Jo 
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From this, it has been shown that the Schwartz inequahty would yield 

< heiO) - 1 < 1, (51) 

and more generally 

\ heir) - I \<\ heiO) - 1 \< I. (52) 

Whenever there is squeezing, these inequalities do not hold for hg{T). Giant violations of 
these inequalities have been predicted for an optical parametric oscillator, and a group of 
atoms in a driven optical cavity, and have been recently observed in the cavity QED system 



Now consider the following quantity 

{lE)^ < {lf{Ef (53) 



After some algebra we find 



< 



< ^^^^HO) (54) 

< 9^'\^)9^^\^) (55) 

In the absence of an external potential g^'^\^) = |C'^f/CgfiP = ^o(O)- ho{0) will be 
nonclassical above 2, we must have bunching to see nonclassical behavior in the conditioned 
fields. Also in the system considered here , we would have h'^ < g^'^'^; when we include an 
optical lattice we have 

^o(O) = (56) 

l^k \^l,g,k\ 

\h.m' = (58) 

I 2^k '^l,g,k\ 



Just looking at the numerator, for two vibronic modes k values, we would violate Eq. (55). 

As with g^'^^r), we obtain an analytic solution using the quantum trajectory method, 
and again we look at weak driving fields. We find 

{a\0)ae{T)a{0)) = (V^, | a, | ^jj,), (59) 
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where | ■0c) is the collapsed state produced by the photodetection event, as in the case of 
g'^'^^r). Once again we need only keep the states with two or less excitations (total in the 
cavity mode or internal energy) for weak driving fields. The result is that 

{n)ss{aeir))cT 



{n)ss{a0{T))ss 
{ae{r))cT 



{aQ{T))ss' 

The expectation value of the field quadrature operator is given by 



(60) 



{ao) = E (V^C;,C„_i,,e-^ + V^CliCn+,,ie^') . (61) 

n,l 

In the weak field limit we have 

(ao) = E (CyCo,^e-^ + C*,Ci,,e^^) . (62) 

I 

So finally then, for weak fields we have 

heir) = ^ °' °' \ (63) 

2^1 {'^1,1 '^0,1 + '^0,1 '^1,1 ) 



For the fluorescent field, we have 



= j^^. (64) 
which in terms of probability amplitudes is 

In Fig. 14 we plot hj^ for g/'-f = 2, /t/7 = 5, for the same choice of four states wc have used. 
We see that in the case of a highly localized atom (equal probability of 20 vibronic levels) 
the nonclassical nature of hj^ is actually enhanced. In the case of an admixture of ground 
and fifth excited states, the behavior of hj^ is relatively unchanged from the ground state 
case. This is due to the insensitivity of hj^ to detunings in the weak coupling limit. In the 
strong coupling regime, as shown in Fig. 15, we see the same general behavior, although for 
the case of 20 equal populations we do sec some dcphasing of the vacuum- Rabi oscillations, 
due to the detunings of the various levels involved. Similar behavior is seen in the case of 
hg^ as shown in Figs. 16 and 17. Note that hg^{0) — 0.0, reflecting the fact that after 
spontaneous emission, the dipole field envelope vanishes. In Fig. 18 we change the level 
spacing. We see that for increasing vibronic level spacing the nature of the nonclassicality 
persists, but there is evidence of the beat frequency between subsequent vibronic levels. 
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BHD signal 

FIG. 13: This is a common experimental setup for measuring hg^r). In this figure, the source 
would be either the transmitted or fluoresced portion of the field. LO denotes Local Oscillator, a 
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T t 

FIG. 14: Plots are /ij^(r) for g/j = 2, k/j = 5 , Ai,+/7 = 0.1. (a) |0) only, (b) ^ [|0) + |5)]. (c) 
Pseudo-Boltzmann. (d) All states equal population. 
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FIG. 15: Plots are /ij^(r) for 5/7 = 3, K/7 = .1 , Ai,+/7 = 0.1. (a) |0) only, (b) ^ [|0) + |5)]. 
(c) Pseudo-Boltzmann. (d) All states equal population. 




FIG. 16: Plots are /iJ^(t) for g/7 = 3, K/7 = .1 , Ai,+/7 = 0.1. (a) |0) only, (b) ^ [|0) + |5)]. 
(c) Pseudo-Boltzmann. (d) All states equal population. 
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FIG. 17: Plots are h^^ir) for 5/7 = 1, k/j = .77 , ^i,+h = O-l- (a) |0) only, (b) ^ [|0) + 
(c) Pseudo-Boltzmann. (d) All states equal population. 




FIG. 18: Plots are /iJ^(t) for 5/7 = 1, K/7 = .77 , Ai,+/7 = 0.1. (a) |0) only, (b) ^ [|0) + |5)]. 
(c) Pseudo-Boltzmann. (d) All states equal population. 
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CONCLUSION 



We have considered the photon statistics of a cavity QED system while including quan- 
tized center of mass motion along the cavity axis. In the limit of weak driving fields we have 
found analytic results for intensity correlations of the transmitted and fluorescent fields; as 
well as for the cross-correlations between the transmitted and fluorescent intensities. We 
find that for intensity correlations for the transmitted field, having a significant population 
outside the ground vibronic level is deleterious to sub-Poissonian statistics, photon anti- 
bunching, and overshoot/undershoot violations. This is explained due to the sensitivity of 
these nonclassical effects to detunings between the atom-cavity system and the driving field. 
It is found that significant population in vibronic levels that are out of resonance by a half 
a linewidth is sufficient to severely modify the results; a highly localized atom, spread over 
many vibronic levels only exhibits nonclassical effects over a very small parameter range. For 
the fluorescent intensity correlations, we do not flnd such a sensitivity, this is due mainly to 
the nature of single atom fluorescence where the atom can only emit one photon at a time. 
The cross-correlations exhibit the assymetry noted by Denisov et. al., and this asymmetry 
is not degraded significantly by a distribution over vibronic levels. 

We have also found analytic results for hg{T) for the transmitted and fiuorescent fields. 
There is no time asymmetry for weak driving fields. The nonclassical behavior in hg{T) is not 
generally degraded by a distribution over vibronic levels; indeed it is sometimes enhanced. 

Future work includes inclusion of 2- and 3-d external trapping potentials, non-harmonic 
potentials, and pressing beyond the weak field limit. 



* Electronic address: 'r icepr @ muohio . edu 
[1] For a comprehensive review, see Optical Coherence and Quantum Optics, L. Mandel and E. 

Wolf, Cambridge (2000). 
[2] E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 (1963). 

[3] Cavity Quantum Electrodynamics, edited by P. Berman, in Advances in Atomic and Molecular 
Physics, Supplement 2, Academic, San Diego, (1994), Any other good reviews that are later??. 

[4] H. J. Carmichael, R. J. Brecha, and P. R. Rice, Optics Communications 82, 73 (1991);R. J. 
Brecha, P. R. Rice, and X. Min, Phys. Rev. A 59, 2392 (1999). 

26 



[5] G. Rempe, R. J. Thompson, R. J. Brecha, W. D. Lee, and H. J. Kimble Phys. Rev. Lett. 67, 
1727 (1991) 

[6] G. T. Foster, S. L. Mielke, and L. A. Orozco Phys. Rev. A 61, 053821 (2000). 
[7] A nice introduction is L Guidoni and P Verkerk, J. Opt. B 1, R23 (1999). 
[8] J. McKeever, A. Boca, A. D. Boozer, J. R. Buck, and H. J. Kimble, Nature (London) 425, 
268 (2003). 

[9] H. J. Carmichael, An Open Systems Approach To Quantum Optics, ( Springer- Verlag, Berlin, 
1993), L. Tian and H. J. Carmichael, Phys. Rev. A. 46, R6801 (1992). 

[10] D. W. Vcrnooy and H. J. Kimble Phys. Rev. A 56, 4287-4295 (1997) 

[11] C. J. Hood and C. Wood, as described by H. J. Kimble et al., in Laser Spectroscopy XIV, 
edited by Rainer Blatt et al. (World Scientific, Singapore, 1999), p. 80. 

[12] H. Katori, T. Ido and M., Kuwata-Gonokami, J. Phys. Soc. Jpn. 68, 2479 (1999), T. , Y. 
Isoya, and H. Katori, Phys. Rev. A 61, 061403 (2000). 

[13] S. J. van Enk, J. McKeever, H. J. Kimble, and J. Ye, Phys. Rev. A 64 013407 (201). 

[14] H. J. Carmichael, H. M.Castro-Beltran, G. T. Foster, and L. A. Orozco, Phys. Rev. Lett. 85, 
1855 (2000). 

[15] G. T. Foster, L. A. Orozco, H. J. Carmichael, and H. M. Castro-Beltran, Phys. Rev. Lett. 85, 

3149 (2000). 



27 



